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Abstract 

Wehrl used Glauber coherent states to define a map from quantum density matrices 
to classical phase space densities and conjectured that for Glauber coherent states the 
mininimum classical entropy would occur for density matrices equal to projectors onto 
coherent states. This was proved by Lieb in 1978 who also extended the conjecture to 
Bloch SU(2) spin-coherent states for every angular momentum J. This conjecture is 
proved here. We also recall our 1991 extension of the Wehrl map to a quantum channel 
from J to K = J + 1/2, J + 1, . . ., with K = oo corresponding to the Wehrl map to 
classical densities. For each J and J < K < oo we show that the minimal output 
entropy for these channels occurs for a J coherent state. We also show that coherent 
states both Glauber and Bloch minimize any concave functional, not just entropy. 



©2012 by the authors. This article may be reproduced in its entirety for non-commercial purposes. 



1 



EHLJPS-September 22, 2012 — coherent states and entropy 



2 



1 Introduction 

Coherent states in the Hilbert space H = L 2 (R n ), which originate in the work of Schrddinger, 
Bargmann, Glauber and others, are certain normalized Gaussian functions parametrized by 
points in classical phase space, (p, g) G K n x W 1 . They are denoted by \p, q). See, e.g., [TJ or 
Section [2] for definitions. 

Given a density matrix p on H (i.e., p is a positive semi-definite operator with trace 
Trp = 1), its von Neumann entropy is S(p) = —Trp In p. This is always non-negative, but 
the usual classical Boltzmann type density (e.g., f(p,q) = Z^ 1 exp — f3(p 2 + V(q)) can have 
an arbitrarily negative entropy — j RnxR n f(p,q) In. f(p,q)dpdq. To rememdy this, and other 
problems with the classical approximation to entropy, A. Wehrl jT3] used coherent states to 
propose another definition, as follows: 

Define p cl {p,q) '■— (p, q\p\p, <?)■ This function has several names. One is the Husimi 
Q-function [5]. Berezin [TJ called it the covariant symbol. The name we shall use here is the 
lower symbol of p (see [6]). Since (27r)~ n J RnxRn \p, q) (p, q\dpdq is the unit operator on H, we 
see that p cl (p,q) is a probability density with respect to the classical measure (2n)~ n dpdq. 
In [U E] it is shown that the trace Tr(/(p)) of a convex function / of p is bounded below 
by the corresponding classical integral (2n)~ n j R „ xRn f(p cl (p,q))dpdq. Together with the 
corresponding upper bound for what is called the upper symbol (or contravariant symbol 
in [1]) they are often referred to as the Berezin-Lieb inequalities. The inequalities are, of 
course, reversed for concave functions. In this paper we shall not be concerned with the 
upper symbol. 

Wehrl uses the lower symbol, to define 



Since < p cl (p, q) < 1, we see that S cl (p) > 0, as desired. 

A question raised by Wehrl is what p minimizes S cl (p), and he conjectured that this occurs 
exactly when p is a one- dimensional projection onto any coherent state, i.e., p = \p, q)(p, q\ 
for any choice of p,q. This conjecture was proved in [7]. Later Carlen [1] gave a different 
proof based on the log-Sobolev inequality; this proof included the uniqueness statement for 
the first time. A decade later Luo [TU] gave a proof based on hypercontractivity, which is 
closely related to the log-Sobolev inequality. 

In [TJ a similar conjecture was made for the coherent states in the Hilbert spaces of the 
irreducible representations of SU(2). These angular momentum coherent states are very 
useful for the physics of quantum spin systems. They are usually called Bloch coherent 




(1) 
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states (21 [TT] and will be the subject of this paper. We shall prove the conjecture in J?J/ 
that the analog of Wehrl's conjecture holds here as well, i.e., that the classical entropy is 
minimized by coherent states. We do not prove that coherent states alone minimize the 
entropy, however. 

Previously, the conjecture had been established only for a few low- dimensional represen- 
tations of SU(2). The 2-dimensional spin 1/2 case is simple as all pure states are coherent 
states. This was already pointed out in [7J. For the 3-dimensional case of spin 1 the conjec- 
ture was solved by Scutaru [13] and by Schupp (12] who also solved it for the 4-dimensional 
representation corresponding to spin 3/2. Bodmann [3] proved a lower bound on the classical 
entropy which is asymptotically correct for large spin J. 

We prove more than that S cl (p) is minimized when p is a projection onto a coherent 
state. We prove this for all concave functions f(t), not just f(t) = — tln(t). In fact, the 
original proof in [7J for the Gaussian Glauber states was for f(t) = —t p , p>l, and the proof 
for — tln(t) followed by taking the limit p — > 1. The extension to f(t) = t p for < p < 1 was 
given by Carlen |4j . To our knowledge it has not been proved for general concave functions. 
In Theorem 12.21 we show exactly that by approximating Glauber coherent states by Bloch 
coherent states in an appropriate limit of large spin. The particulars of this approximation 
are in the appendix. 

In order to prove the conjecture for all spin we utilize a generalization of coherent states, 
called coherent operators introduced by us in [8]. These are operators that map density 
matrices in one SU(2) space, characterized by an angular momentum (or "spin") J to a 
density matrix in a spin K space. These maps are, in fact, quantum channels, i.e., completely 
positive trace preserving maps, as will be made clear in Lemma 14.11 

The Bloch coherent states map density matrices from J to functions on the classical 
phase space, i.e., the 2-sphere § 2 (which can be thought of as K = oo [5J) via the lower 
symbol map mentioned above. We do this in small steps, so to speak, by going from J to 
J + 1/2 to J + 1, . . .. For each finite K on the way we prove that projections onto coherent 
states in J minimize the von Neumann entropy of the lifted density matrix in K. In other 
words we determine the minimal output entropy of the quantum coherent operator channels 
to be the entropy of the output of a coherent state. We then show that after an appropriate 
scaling, the limit K — > oo gives us the desired classical (lower symbol) entropy, and thus we 
prove the conjecture for the entropy and for any concave function. 

An important observation in our procedure is to note that the quantum coherent operator 
channels have a simple expression in terms of bosonic second quantization, i.e., bosonic 
creation and annihilation operators. 
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Coherent states can be generalized to any compact semisimple Lie group, not just SU(2) 
(see pj)J EE]), and there we expect that a similar result holds. 

In Section [2] we define Bloch coherent states, lower symbols, and discuss the corresponding 
Berezin-Lieb inequality. We also introduce the quantum coherent operator channels from 
J to K. (While we are interested here in K > J, the map is also defined when K < J.) 
In Section [3] we derive a more explicit formula for the quantum coherent operator channel 
which allows us, in Section HI to give a bosonic second quantization representation of the 
channels. In fact, Section [3] is not important for our main conclusion as we could have 
defined the quantum channels from the second quantization formulation in Section HJ We 
include Section [3] in order to connect to our previous work in [8J. Following that, we show, 
in Section [5j that coherent states in J minimize the output von Neumann entropy in K or 
more generally the trace of any concave function. In the last Section [6] we study the classical 
limit K — > oo and use the Berezin-Lieb inequality to prove the conjecture on the classical 
entropy. 

Acknowledgment: Part of this work was carried out at the Mathematics Department 
of the Technical University, Berlin and at the Newton Institute, Cambridge. We are grateful 
to both and, in particular, to Ruedi Seiler for hosting our stay in Berlin. 

Thanks go to Eric Carlen, Rupert Frank, and Peter Schupp for their helpful comments 
on a preliminary version of the manuscript and to Anna Vershynina for a careful reading of 
the manuscript. 

This work was partially supported by U.S. National Science Foundation grant PHY- 
0965859 (Elliott Lieb ), by a grant from the Danish research council (Jan Philip Solovej) 
and by a grant from the Simons Foundation (#230207 to Elliott Lieb). 

2 Basic definitions and main results 

For all integer or half- integer J we let %j denote the spin J representation space of SU(2), 
i.e., %j = C 2J+1 . The corresponding classical phase space is § 2 , the unit sphere in M. 3 . For 
each point u G § 2 we have the one-dimensional coherent state projection = j(oj\ 
projecting Tij onto the subspace of maximal spin in the direction u, i.e., the one-dimesional 
subspace of T-Lj corresponding to the eigenspace of u ■ Sj with eigenvalue J. Here Sj is 
the vector of spin operators, i.e, the representation on T-tj of the standard generators S = 
(S x , S y , Sg) of SU(2). The vector \co)j is only defined up to a phase, but this will not play a 
role here as only the projection is important. We will use the notation that 7, S 2 are 
respectively the north and south pole. 
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The coherent state transform is based on the identity that 
2t/ + 1 f , . , , , 2J+1 [ . 



where Ij is the identity on %j. If p is a density matrix on Hj its lower symbol is the function 
on § 2 given by 

$°°(p)M = >|p|u;) J = Tr J (Pjp), (3) 
where Trj is the trace on Hj. The classical entropy of p is 

S C \ P ) = j— / $°°(p)(u)H$°°(p)(u))du. 

We are using the notation $°° for the lower symbol since we shall consider it as the natural 
classical limit k — > oo of the quantum channels $ fc defined below. 

The Berezin-Lieb [HE] inequality for the lower symbol states that for any concave function 
/ : [0, 1] -»• R we have 

Trjf(p) < [ f{*°°(p){ U ))du,. (4) 



47T 

The inequality follows from (j2J) as a consequence of Jensen's inequality. 

The conjecture from j7] that we shall prove here is that 5* 01 is minimized when the density 
matrix is any coherent state projection, e.g., p = | T)j/(T I- I n this case the lower symbol is 
$°°(| t)jj(t |)( w ) = T)j| 2 - m fact, we shall prove the more general statement that the 
same is true if the function — tln(t) is replaced by any concave funtion. Our main theorem 
is the following. 

2.1 THEOREM (Lower symbols of Bloch coherent states minimize concave averages). Let 
f : [0, 1] — > R be a concave function^. Then for any density matrix p on Tij we have 

2 7-4-1 /* 2 7-1-1 f 

' fWpW*" > — A i- / /(|>|t>/)^. (5) 



47T ./s2 ~ 47T 

By SU(2) invariance t could be replaced by any other point on S 2 . 

The following analogous result for the Glauber coherent states is proved by an easy 
limiting argument which we give in the appendix. 

2.2 THEOREM (Lower symbols of Glauber coherent states minimize concave averages). 
Let f : [0, 1] — > R be a continuous concave function with /(0) = 0. Then for any density 



1 It is, in fact, enough to assume that / : [0, 1) — > R, i.e., to allow that lim t ^o- f{t) — — oo. Only coherent 
state projections have lower symbols that attain the value 1. If p is not a coherent state projection we can 
find a concave function / : [0, 1] — > R such that / > / and f { j{uj\p\uj)j) — f{ J {uj\p\uj) J ). 
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matrix p on L 2 (M. n ) we have 

(2ir)~ n / f( P cl (p,q))dpdq > (2tt)-» f f(\( Po ,qo\p,q}\ 2 )dpdq (6) 
for all po,qo G M n . 

2.3 REMARK. As in the main Theorem 12.11 we could have allowed lim^i- f(t) = — oo (see 
Footnoted]). We could, in fact, also allow /(0) 7^ 0, but in this case the integrals on both 
sides of ([6]) are either both +00 or both —00. Even if /(0) = the integrals may still be 
+00, but the inequality holds in the sense that either both sides are +00 or the right side is 
finite. 

Using the fact that the Glauber coherent state \p,q) E L 2 {R n ) is explicitly given by 

7r~ n/4 exp(-(x - q) 2 /2 + ipx), 

we have 

(p, q\ip) = 7r~ n / 4 J exp(— {x — q) 2 /2 — ipx)ip(x)dx, 
for ip G L 2 (R™). The inequality ([6]) for the rank one state then states that 

/ f(\(p,m\ 2 )dpdq 

is minimized for concave / when ip is a Glauber coherent state. 

We now define the quantum coherent operator channels. We refer to [H] for details. For 
fixed K and J we let P_ be the projection in %j <g> Hk onto the minimal total spin \K — J\, 
i.e., onto the unique copy of U^ A QHjVHk on which the tensor product representation 
acts irreducibljo For simplicity we omit in our notation the dependence of P_ on K and J. 

In the language of elementary quantum mechanics a particle of angular momentum K 
and one of angular momentum J can combine in exactly one way to produce a composite 
particle of angular momentum \K — J\. The Hilbert space of this composite particle is the 
subspace %\k-j\ Q Hj ® Hr- 

If we let k = IK — 2 J 6 Z we consider the map $ fc from operators on %j to operators 
on %k defined by the partial trace 

97x1 

®\p) = 2 \ K _ j\ + 1 TTj{P - {p ® Ik)) - (7) 
2 Strictly speaking the isometric imbedding of H\k-j\ m to H.j®1-Lk is given uniquely only up to a phase. 
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This is a trace preserving completely positive map (see also ( TT6l) and Lemma T4. II) . i.e., using 
the language of quantum information theory it is a quantum channel. The trace preserving 
property is easily seen since the partial traces Tr jP_ and Tr^-P- are both proportional to the 
identities. In particular, $ fc maps density matrices to density matrices. In the notation we 
have for simplicity omitted the dependence of $ fc on K and J and only kept the dependence 
on the difference in dimension k = 2K — 2 J. 

Our main result about these channels is that they are majorized by coherent states in 
the following sense. 

2.4 THEOREM (Coherent states majorize $ fe ). For a density matrix p on Hj and k = 
2(K — J) the sequence of eigenvalues of the density matrix $ fc (p) is majorized by the sequence 
of eigenvalues of § k (\u) (u)\), which by SU(2) invariance is independent of u G § 2 . 

To say that a finite real sequence a\ > 02 > • • • > a a/ majorizes another real sequence 
b\ > b 2 > ■ ■ • > b M , written (a 1; . . . , clm) >~ (&i, • • • , &m) means that 

a\ >b\ 
a x + a 2 >b\ + b 2 

a x + . . . + a M -i >bx + ... + b M -i 
ax + ■ ■ ■ + a M —bx + ■ ■ ■ + b M - (8) 

Note the equality in the last condition (jSJ). 

It is a fact that (ax, ■ ■ ■ , ajw) >~ (bx, ■ ■ ■ , &m) if and only if 

M M 

3=1 3=1 

for all concave functions / : R — >■ R. This is often called Karamata's Theoremjf] cf. [H] 
Remark 4.7 after eq. (4.5.4). It is, in fact, enough that the concave function is defined on 
the interval [am, ai]- 

If A and B are two Hermitian matrices of the same size we write A y B if the eigen- 
value sequence of A majorizes the eigenvalue sequence of B. The notion of majorization of 
sequences can be easily generalized to infinite summable sequences and trace class operators 
but we will not need this here. 

3 In [9] the concave function is assumed to be monotone increasing. With the assumption of equality in 
([5]) the assumption of monotonicity is not required. 
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It is an easy exercise, using the variational principle, to prove that if A, B, C are Hermitian 
matrices such that A y B and Ay C then 

A y XB + (1 - X)C, (9) 

for all < A < 1. 

As a consequence it follows from Theorem 12.41 that the minimal output (von Neumann) 
entropy of the channel $ fc is achieved for a coherent state, i.e., 

minS($ fe (p)) = S($ k (\tu)(tu\)). 
p 

More generally, the output of coherent states minmize the trace of concave functions. 

2.5 COROLLARY (Minimization of the trace of concave functions). If f : [0,1] — > M is 
concave and p is a density matrix on Hj and k = 2(K — J) then 

Tr K [f(^p))]>Tr K [f(^(\u)(u\))] 

for all w£§ 2 . 

Of course, the inequalities about concave functions are reversed for convex functions. 

3 A formula for P_ 

Our goal here is to find an explicit formula for the projection P_ that projects H j £g> H-k 
onto the subspace Hk-j, under the assumption that K > J. 
We start by choosing the standard preferred basis 

\M) L , M = -L,...,L, 

in Hi, in which S z is diagonal and S x is real. This specifies the basis up to an over-all phase. 
We introduce the anii-unitary map Ul '■ Hi —> He given by 

L L 

u L oc m \m) l = {-i) L - M Tm\-M) L . 

M=-L M=-L 



This map has the property that 

U L l S L U L = -S L (10) 
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It is the unitary e lnSy followed by complex conjugation in the preferred basis. Any anti- 
unitary satisfying ( flOj) agrees with Ul up to an over-all phase. Note that U L l = (—1) 2L Ul 
and hence 

(u L m = {u L ^\u L u L( t>) = (-i) 2L (u L m, (ii) 

for all 0, if) G Hl- 

If K > J then %k Q Hk-j <8> %j (recall that Wk-j Q Hj <8> %k as we said in the 
beginning). We thus have a sesqui-linear map 

where the partial inner product j^ll^L is defined by the inner product in Hk-j as follows 

k-j(v\j(^U)k)k-j = { k- J)9 j( t 1®iI>\<I>\k-j) 9 j 
for all r] G Hk-j, where the last inner product is in Hk-j ® 

3.1 THEOREM (Formula for P_). If K > J then we have for all ip G Hj and G H K 

P-(4>®<p) = t i JjJMfa, (12) 

where fi G C satisfies 

_ 2(g - J) + 1 
2fC+l 

Proof. Formula f)12p above for P-(if) <8> 0) is clearly a bilinear map in -0 and 0. It is thus 
enough to prove the above formula for a linear spanning set for and ip. Such spanning sets 
are provided by the coherent states and \uj) k for uj G S 2 . It is thus enough to prove the 
formula for ip = \uj')j and = For simplicity we will write Uj\uj') t = \Ujuj'\ 7 , where Uj 
is the anti-unitary defined above. 

In the following we let Sj and Sk, respectively, denote the spin operators on %j and Hk, 
respectively, and we let S be the total spin operator on 'Hj®'Hk, i-e., S = Sj®Ik + Ij®Sk- 
Since uj ■ Sk\w) k = K\uj) k , and since rj = \uj')j — \U jcj) jj (U i juj\uj')j has no component in the 
subspace uj ■ Sj = — J, it is clear that P_r\ ® \uj) k = 0, for otherwise the total uj ■ S component 
of this vector would be bigger than the maximal possible namely K — J. Hence 

P^co'l ® \uj) k = JJJM<S)jP-\UMj ® \")k- (is) 
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We have uj ■ S\Uju)j <g> \u) K — (K — J)\Ujuj)j (g> and thus, since P_ commutes with uj ■ S, 

P-IZ/ja;), <8> \u) K = fi! \uj) k _ 3 = n' Jp\\u) K (14) 

for some complex scalar //. Here we have used that j(uj\\uj) k = \oS) K _j in %k-j ® Uj, since 
\oj) k = \uj) k _j <8> Inserting (TH|) into (IT5|) we obtain 

Since Z7J = (— 1) 2J and Uj is anti-unitary we have Uj(\u')j — rj) — {—1) 2J j {oj'\Ujuj) j \uj) j . 
Moreover, Ujt] has no component in the space uj ■ Sj = J, hence j(C^j?7||w) k = and thus 

with /i = (— l) 2J yu', which is what we wanted to prove. 

We can find the modulus of \i from the fact that with k = 2(K — J) is trace 

preserving 

M 2 =(.,<C^I ® >|)P-(|!7^ ® k)J = Tl V (lV, { (P-k)„->l)) 

- 2 ( if - J ) + 1 W M,>|)- 2(Jr - J) + 1 



□ 



If k = 2(K — J) > we therefore have 



2 7+1 

.(01^(1^^1)1^ = 2^iIUWII^IIk-j- (is) 

If we introduce the channel 

&{p) = ^(UjpUj 1 ) (16) 

we see that 

~ 2 7+1 

or equivalently 

~ 2 7+1 

^{p) = -^ 1 Pk(Ik-j9p)Pk 1 (17) 
where Pk is the projection onto %k in Hk-j ®7~Lj- In particular, $° is the identity map. 
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If K < J the corresponding result is that 

p^®4> = n' K (u K <t>Ul, \^ 2 = 2{J 2j+\ +1 - (18) 

In particular, in this case 

M^ k Km,Am<p) K = uuk4>ww 2 j-k. (19) 

Hence if K < J and we now set <lH fc l(p) = U K 1 ^ k (p)U K we obtain 



K m~ lkl mM\)\<p) K = \\ K 



4 The Bosonic formulation 



r\U-K- 



The space Hj may be identified with the completely symmetric subspace 0gy m "Hi/ 2 of the 
tensor product ® 2J "%i/2- 

A particularly simple way to see this is to use the Schwinger representation of spin 
operators in terms of creation and annihilation operators. Let %\/2 be the one-particle space 
and let and aj be the creation operators corresponding to spin up and down respectively. 
They are the operators which, for all positive integers £, map <S>f ym ^i/2 to ®fym^i/2> sucn 
that for Ve<8>f ym ?*i/2 

a *V = VI + IPsymd t)i ® 

and likewise for a|, where P sym is the projection onto the symmetric space ®sym^V 2 - The 
annihilation operators a^, a± are the adjoints of a| and a|. 

The symmetric subspace of (£) 2J W.1/2 is the subspace corresponding to 2 J bosonic parti- 
cles, i.e., the subspace 

a^a^ + a|a^ = 2 J. 
On this 2 J + 1-dimensional subspace we observe that the operators 

S x = ^(a^ + a|a t ) 
5 W = ^-Xa\a x - a* x a t ) 
S z = ^(a* t a t - aja 4 ) 

satisfy the correct commutation relations and S% + Sy + S% = J(J + 1). 
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The spin representation on T-Lj may then be identified with the space of 2 J bosons over 
a 2-dimensional one-particle space. In particular, the coherent state G TLj is in the 
bosonic language the pure condensate wave function ((2 J)!)" 1/2 (a* ) 2J |0), where |0) is the 
vacuum state (i.e., the state of zero particles) and a* is the creation of a particle in the state 
\u>)±, i.e., a* = {\ \cu\al + ,(I \uj\a1. We will use the canonical commutation relations 

2 2 2 1 2 2 + 

that all creation operators commute and [a<y,a*] = 1 (uj'\uj) 1 . This gives, in particular, 

2 2 

a^a| + aj,a| = a^a-j- + a^a^ + 2 = 2 J + 2 (20) 

on "Hj. 

The channel $ fc has a simple form using creation and annihilation operators. 

4.1 LEMMA (The channel $ fc in second quantization). If p is a density matrix on Hj = 
® sym "^1/2 an d K = J + I / or some integer k > i/ien 

$fc(p) = ^TT)T <---<P Q n---^- ( 21 ) 

1 J ' ii,...,i fe =U 

// K = J — | wi/i Ic > we /icwe 

^" (p) = i7Tr a ik ---a il pa* 1 ---a* k . (22) 

Proof. We use the expression ([17]) for The space = (S^yL^i/ 2 * s the totally sym- 
metric subspace of (^i^T^ H1/2 J ® (®sym ^1/2) • Thus by the definition of the creation 
and annihilation operators^ 

~ 2J + 1 

$ fc (p) = , ■ : V ® ®p® x<*l| ® •••® (23) 

2il + 1 Z ' 2 2 2 2 



«iv4=t,J- 
2J+1 (2 J) 



oJ Jk ---o* 1 po<i---Oi k - (24) 



2K+1 (2101 

v ; h,...i k =t,l 

The case K < J follows in the same way. □ 

4-2 REMARK. Note that ( }2~li) and ( 122|) are the Kraus representations of the completely 
positive trace preserving maps 

We shall use this to calculate the action of $ fc on coherent states. 



4 The meaning of the operator <S> p® (ip\ in (P2IJ1) is clear if p is a rank one projection and it is 

defined in general by linearity 
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4.3 LEMMA (The action of $ fc on coherent states). If K = J + | /or some integer k > 
i/ien $ fc (| t)jj(t I) ^ as ^ e orthonormalized eigenf unctions 

<f* = (j!)- 1/2 ((2 J + fc - j)!)" 1/2 («I)'(4) 2 j = 0, . . . , 2 J + A; = 2A" 

ii/ii/i corresponding eigenvalues 

c>k 2J + 1 k\(2J + k-j)l 2J + 1 (2(K - J))\(2K - j)\ 
3 2J + k + l{2J + k)l{k-j)\ 2K + 1 (2K)\((2(K — J) — j)V 

for j — 0, . . . , k and zero for j > k. Note that the eigenvalues are listed in decreasing order. 

Proof. This follows immediately from Lemma 14.11 since 

$ fc (l tL(t I) - (2 j+t+ 1)! ]C Q)((2^)0- 1 (4) i (4) 2J+fc 'i Q )(Q]4 WQ l- 

□ 

An important ingredient in the proof of our main Theorem 12.41 below is to study the 
operator 

m 

T-iC,fc+l I iC,k+l\ /iC,k+l\ * . I ,C,k+l\ I /C,fe+l| * 

j=0 

for m < k. Using the fact that a|0^' fc+1 = \fj (jyf^i for j > 1 and a^(j)^ k+1 = and 
a t^j = + k + 1 — j (j)f' k for j = 0, . . . , m we find that 

m— 1 

r£ fc+1 = $>J + fc + 2)\4>f h )(<j>f k \ + (2J + k + l- m)\^ k )(^ k \. (25) 

3=0 

4-4 REMARK. The expression in Lemma 14.11 for the channel may be generalized to 
define analogous channels between bosonic many-particle spaces where the one-particle space 
instead of being 2-dimensional, as T-Li/2, could be of arbitrary finite dimension. We conjecture 
that Theorem 12.41 holds also in this case in the sense that pure condensates majorize these 
channels. 



5 Proof of the main theorem for the channels <E> 



k 



Proof of Theorem\2.4\ If k — 2(K — J) < then $ (|o;)(a;|) is rank one and the result is 



obvious. We now consider the case k = 2{K — J) > 0. We first point out that from (J2J) it 
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is enough to consider the rank one case, i.e., p = \if>)(if>\ for some if) G Hj. Since Ir-j has 
rank 2{K — J) + 1 = k + 1 it is clear from f fTTl) that <3> fc (\if)){if> |) has rank at most + 1. 
It is of course equivalent to consider the channel <3> fc . Let \j(if>), j = 0, 1, . . . , 2J + k be 
the eigenvalues of $ (IV'KV'I) i n decreasing order and counted with multiplicity. Let <fq, 
j = 0, 1, . . . , 2 J + k be the corresponding orthonormalized eigenvectors. For all m > we 
have ^jLo^jK^) = Tr = 1. The claim is that, moreover, 

m m 

£^)<£Af, (26) 

3=0 3=0 

for m — 0, . . . , A; — 1, where A^ ,A; are the eigenvalues for the coherent states, which were given 
in Lemma [4. 31 

We shall prove ( |26l) by induction on m. For m = this is easy since we clearly have from 



( ITTjl that Aq^) < 2^+1 = ^o'^- us now assume that we have proved fl2"B"|) for all integers 
up to m — 1 for some m > 1 . We want to prove it for m. We shall do this by induction on 
k. For k < m, (1261) is an equality since both sides are 1. Let us assume, therefore, that we 
have proved (|26|) up to some k > m. We want to prove it for k + 1. 

Since 0o +1 , . . . , G Hj + (fc + i)/2 are the orthonormal eigenvectors corresponding to the 
m top eigenvalues of $ fe+1 (|?/>) (if>\) we have 

m for , i\| m 

i=o 1 j=o ti,...,i fc+1 =t4 

(2J+1)! 
(2J+A; + 2)!. 

1 1 2J+fc 



where 



r = E Kl^fWK + <hI# +1 >(# +1 K) , 

3=0 

is an operator on the space Hj+k/2- Observe that since are 2 J + k + 1 particle states 
we have 

m 

Trr = E(0i +1 K«t + oKI^ 1 ) = ( m + 1)(2J + k + 1). 

3=0 
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J+fc/2" 



Likewise, we have from ( 1201) with 2 J replaced by 2 J + k, the operator inequalities 

m m 

< r =a r J2 + a+E < «t4 + = (2 J + A; + 2)/ w 

i=o i=o 

We would therefore get an upper bound to the expression in ( 1271) if T is replaced by 

m— 1 

£(2 J + k + 2)\^)(4> k j \ + (2J + Ac + 1 - m)|0* W* |. (2* 
i=o 

This gives the bound 

i=o i=o 

2J + Hl-mA , m+1 



2J + /C + 2 f-< jvr/ 2J + /C + 2 . 

We conclude from the induction hypotheses on both m and k that 

m o 7 i 7 i 1 m l 1 m_1 m 

j w ^ 27+^+2 E A , + 2J+fc+2 E A , = E A i 

J'=0 i=0 j=0 j=0 

That the last recursive identity holds for the coherent eigenvalues follows since for coherent 
states we know from (1250 that F is, in fact, equal to the optimizing expression (12 8j) . It can 
also be seen from the explicit formulas in Lemma [4.31 (e.g. using induction). The induction 
is thus complete. □ 

Note that as a special case we have seen in the proof that for k = 2(K — J) > 

2 7-1-1 

< ^ (29) 

6 The classical limit of the channels $> k 

In this section we consider the limit as the dimension of the output space tends to infinity, 
i.e., K — > oo. 

6.1 LEMMA (The large K limit of coherent state outputs). Let f : [0, 1] -+ R be a 
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continuous function. Then 



2J+1 



lim 

K^oo 2K + 1 



TV 



K 



f 



2K + 1 
2J+ 1 



$ fc (l TUT I) 



2J 



47T 



/(|>|UI 2 )^. 



(30) 



Proof. We may of course replace the channel $ fc by $ fc . This a simple calculation based 
on the explicit expressions Lemma [4.31 If 9 denotes the polar angle of u, i.e., cos(#) is the 
^-component of u then \j(u)\ t}/| 2 = cos 4 J (0/2). The integral over the sphere may hence be 
rewritten as 



2J + 1 

47T 



/(|>|t) 7 | 2 )^ = (2J+l) / f{t 2J )dt. 

'o 



On the other hand the explicit eigenvalues in Lemma 14.31 give 



2J+ 1 
2K+1 



Tr, 



./ 



2K + 1 
2J + 



r **(i tur I) 



3=0 V V 



2K + 



27 1 2{ ^ J) n f (2(K - J))\(2K - j)\ 



2K + 



3=0 V 



(2K)\((2(K- J)-j)\ 



It is an easy exercise which we leave to the reader to show that this converges to the above 
integral in the limit as K — > oo. □ 

It can be proved that the same limiting equality ( 1301) holds even if | t)jj(T I is replaced by 
any density matrix p on Hj, at least for a large class of functions /. We will not do this here. 
Instead, we shall restrict ourselves to an inequality similar to (130]) for concave functions /. 
To do this we shall use the Berezin-Lieb inequality (jl]). 

6.2 LEMMA (The classical integral dominate the trace of concave functions of $). Assume 
f : [0, 1] — > R is a concave function. Then for any density matrix p on T-Lj we have for all 
integers k = 2{K — J) > 



2J + 1 
2K+ 1 



TV 



K 



f 



{2K + 1,. . 



V 2J 



+ 1 



- 2J i + 1 / f(Mp\ u )j) du - 

4vr J§2 



Proof. Again we consider the equivalent channel The result follows from the Berezin-Lieb 
inequality 01]) if we can show that the lower symbol of $ fc (p) satisfies 



M<b k {p)\u) K 

This is straightforward from ffTT]) . 



2J + 1 
2K + 1 



j{u\p\u)j- 



□ 
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We are now in a position to prove the main result Theorem 12.11 In fact it is the analog 
of the main Theorem 12.41 or rather the equivalent formulation Corollary 12.51 for the classical 
map $°° from density matrices on %j to functions on (the classical phase space) S> 2 . For 
classical functions the trace is replaced with the integral over phase space. 



Proof of Theorem \2.1[ From Corollary 12.51 and Lemma 16.21 we have for all integers k 
2{K - J) > 



2J + 1 
2K + 1 



Tr 



K 



f ( 2K 



< 



2J+1 
2J+ 1 

4:71 



■**(ltUM) 



< 



2J + 1 
2K + 1 



Tr 



K 



f 



{2K + 1,. . 

$ k (p) 



2J 



f(j(u\p\w)j)du- 



The result now follows from Lemma 16.11 



□ 



A Proof of the generalized Wehrl conjecture, Theo- 
rem 12.2 



It is enough to prove Theorem 12.21 for n = 1. The general case follows by induction as follows. 
Assume we have proved it for n— 1. Then for each (p', q') e R n_1 we define an operator py,?' 
on L 2 (R) by 

WttrfW) = (01 ® (p',q'\p\p',q') ® Wi- 

Then 

Pp',q' = (Trx 2 (K)Pp',?') Pp',g' 

is a density matrix on L 2 (IR) and we get from the inequality for n = 1 that 

l X 



2<T" / /(p cl (p,g))^g>(27r)-" / f(Tr Lm p p , >q ,\(0,0\ Pn ,q n )\*)dpdq. 



We have, however, that Tr^^joy^/ = (p', g'lTr^pjp', g'), where Tr n p is the density matrix on 
L 2 (IR n_1 ) obtained by taking the partial trace on the n-th variable. Thus, from the induction 
hypothesis, we find 

/ f(p d (p,q))dpdq >(27c)- n [ f M/n n,../ ./M2,/n m.. " 

>xR n Vk«x1 



(2tt)-" / /(p cl (p,g))«>(27r)-" / /(|(0,0b',g')rl(0,0b„,g n )|^)^g 

<27r)- n / /(|(0,0|p,g>| 2 )dpdg. 



"yD" 



It remains to prove §6§ in the case n = 1. We first observe that, by possibly replacing 
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f(t) by f(t) + at, we can assume that / is non-negative. Moreover, using the monotone 
convergence Theorem we can assume that / is piecewise linear. In this case we have the 
inequality 

!/(*)- /(y)|<C7|x-y| 
for some C > and all x,y G [0, 1]. Hence for all density matrices pi, pi we have 

\M(P, ?)) - f(pf(P, Q))\dpdq <C [ \pf(p, q) - pf(j>, q)\dpdq 

l /K"xR n 

<C (P,q\\Pi~ P2\\p,q)dpdq 

Jm. n xR n 

=C\\pi - p 2 1| i, 

where the the norm on the right is the trace norm. Hence it is enough to prove ([6]) for a 
subset of density matrices that is dense in trace norm. 

Let |n) denote the eigenfunctions of the harmonic oscillator 

d 2 \ 

+ q 2 | n ) = (2n + l)|n). 



We will prove ([6]) for the dense family of p satisfying the property that there exists a positive 
integer N such that 

(n\p\m) = if n > N or m > N. (31) 

We shall apply the convenient complex notation, where z = 2^ 1 ^ 2 (q+ip) and ~z = 2~ 1 ^ 2 (q—ip). 
The Glauber coherent states may then be written 

2 io Z 

n). 



\p,q) = \z) = J2e-^ 2 ^= 

n=0 



We identify the subspace Hj with span{|0), . . . , |2J)} in L 2 (R) in such a way that \M) 3 = 
\M + J) for M = — J, . . . , J. Moreover, we also identify the 3-sphere with the complex plane 
through stereographic projection, such that the measure is (1 + \z\ 2 / 4)~ 2 d 2 z. With these 
identification we can conveniently write the Bloch coherent states (see [B]) as 

2J /0 T \ 1/2 



(1 T\ 1 

}, = E J (l + \z\ 2 /4r J (z/2) n \n)- 
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It is now straightforward to see that if p satisfies assumption ( 13TI) then 

(z\p\z) 

Using the fact that (1 + \z\ 2 /(2J)) 
easily see that for J > N + 2 

,<(2/J)*Z| p |(2/ J)^ < C^(l + M7(2(iV + 2)))- 2 

if p satisfies (|3"T1) . Since / is non-negative and bounded above by t h-> at for some a > we 
immediately find from dominated convergence that 

/ f((z\p\z))d 2 z = lim / /G((2/J)^|p|(2/J)^) J )(l + |z| 2 /(2J))- 2 ^ 

= lim ^ / /(>|p|4)(l + M74)- 2 d 2 z. (32) 
Our main Theorem 12.11 and the observation that |0) = \0)j implies that 

/ /G(z|p|4)(l + |z| 2 /4)- 2 ^> / /(| J (z|0)| 2 )(l+|^| 2 /4)-Vz. 

Since the density matrix |0)(0| clearly satisfies (13~T1) we see from our main result and ( 1321) 
that ([6]) holds for all p satisfying fl3T|) and hence by approximation for all density matrices 
on L 2 (R). □ 
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